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ABSTRACT 

Wrapping around the Milky Way, the Sagittarius stream is the dominant substructure in the halo. 
Our statistical selection method has allowed us to identify 106 highly likely members of the Sagittarius 
stream. Spectroscopic analysis of metallicity and kinematics of all members provides us with a new 
mapping of the Sagittarius stream. We find correspondence between the velocity distribution of stream 
stars and those computed for a triaxial model of the Milky Way dark matter halo. The Sagittarius 
trailing arm exhibits a metallicity gradient, ranging from —0.59 dex to —0.97 dex over 142°. This 
is consistent with the scenario of tidal disruption from a progenitor dwarf galaxy that possessed an 
internal metallicity gradient. We note high metallicity dispersion in the leading arm, causing a lack 
of detectable gradient and possibly indicating orbital phase mixing. We additionally report on a 
potential detection of the Sextans dwarf spheroidal in our data. 
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Subject headings: galaxies: dwarf, Galaxy: formation, Galaxy: halo, Galaxy: kinematics and dynamics 


1. INTRODUCTION 


The Sagit tarius (Sgr) dwarf was discovered by Ibata 


et al. (1994) and its associated stream was mapped us¬ 
ing M-gi ant stars in the Tw o Micron All Sky Survey 


(2MA SS;[Skrutskie et al. 2006|) catalog by Majewski et al. 
(2003). T'his stream traces an arc which wraps at least 


twice around our Milky Way galaxy, passing through the 
gravitational potential of the Milky Way. This stream 
can h elp us probe the da rk matter halo around the Milky 


Way (Ibata et al. 2001), as well as understand the pro- 
genitor of the current dwarf and stream system. In par¬ 
ticular, the leading arm of Sgr is known to be sensitive 


matter halo (e.g., 

Ibata et al.|2002 Majewski et al.|2004 

Vivas et al.||2005 

Law et al.|[2U09p. 


by Gasey et al. (2012aI. They found th at the Milky 


Way t r iaxial da.rk matter halo model of Law fc Ma- 
jewski (2010) allowed them to improve their match be¬ 


tween modeling and data for the Sgr stream. The tri¬ 
axial model has since gained further support with the 
data from blue hor izontal branch stars in the Sgr stream 
(Yanny et al.|2009), where velocities closely match those 
predicted in the triaxial model case. Although several at- 
tempts have been made to map the stream (for example, 
Majewski et al.||2003[ [Johnston et al.||2012] [Belokuro^ 
et al. 2014 Pila-L )iez et al.||2014p there are still conflicts 


on M nFy Way halo shapes determined from Sgr models 

Johnston et al.|2005||Law & Majewski 

- ^ -- 


e.g., Helmi|2004 [ 

i'issera et al. 


Mol 


eg fe Widrow|201^' Vera-G^iro 


.mi ||20l3|1TEata et al.||2013p). 

Une particular feature that models have had some issue 
reproducing i n the Sgr stream is an ap parent bifurcation 
discovered by Belokurov et al. (2006). This bifurcated 
nature of the str eam has been conhrined in both Galac- 
tic hemispheres (Koposov et al.||2012 Slater et al.||2012| 
and was pro posed to result I r oin a spherical dark mat- 
ter halo by Fellhauer e t al. (2006). It was shown by 
Peharrubia et al. 120I0|) that a bifurcation can naturally 
arrive from mode l s usin g a disky progenitor for Sgr, but 
Peharrubia et al. (2011) found none of the intern al rota¬ 
tion that wo uldoe indicative of such a scenario. |Newby| 
et al. (2013) have suggested that the two tails of the Sgr 
stream may have diffe rent progenit o rs. An additional 
possibility explored by Casey et al. (2012a) is that the 
bifurcation may not need to be ht to current models if it 
was caused by some dramatic recent event. H owever, re¬ 
cent tracing of several Sgr stream members by [Belokurov 


et al. (2014) led them to conclude that no existing simu- 
lation tor disruption could explain the observational data 
and the problem has yet to be solved . 


Simulations by Bailin & Steinmetz (|2005|) prefer triax¬ 
ial dark matter haloes for galaxies when producing large 
scale structure formations but it may not necessarily be 


E.Hyde@uws.edu.au 

^ Macquarie University, Physics and Astronomy, NSW 2109, 
Australia 

^ Australian Astronomical Observatory, PO Box 296 Epping, 
NSW 1710, Australia 


required in the case of the Sgr dwarf ( Ibata et al.||2013 ). 
Although there is some evidence for a Milky Way triaxial 
dark matter halo model, there are still numerous ques¬ 
tions to be answered about the progenitor of the stream, 
in particular its metallicity. Previous mapping of the 
core and stream of Sgr have included estimates of the 
metallicity distribution function (MDF). However, areas 
which show a mixed orbital phase and the gradients in 
metallicity with distance along the stream remain uncon¬ 
vincingly mapped. The mixed orbital phase regions oc¬ 
cur when debris or stars stripped from the satellite orbit 
land within a range of azimuthal time periods about the 
satellite’s own. This scenario would lead to phase mixing 
of d ebris ahea d and behind the satellite along its orbit, 
and Johnston (1998) predicts that the tidal streamers 
from the Sagittarius dwarf galaxy currently extend over 
more than 27r in azimuth along their orbits. The range 
extended by the debris at different points would increase 
the type or age of stars present at a given point in the 
stream. For this reason some areas of the stream are pre¬ 
dicted to have a wider range in metallicity values than 
others ([Law et al.||2005|). 

The population near the core of Sgr has be en found 
to have a mean value of [Fe/H] ^ —0.5 dex (lyacciar i 


et al. 2002 Bonifacio et al. 2004 Monaco et aT. 2005). 


i'his value is known to vary as we move out from the 
core, potentially reflecting the properties of the original 
progenitor of the dwarf and stream. 

In the leading arm of the Sgr stream, jChou et al.|(|2007|) 


found a mean metallicity of ([Fe/H]) = —0.72 dex. Gasey 


et al. (2012a I find K-giant members of Sgr which are no- 
tably more metal-poor (([Fe/H]) = —1.7±0.3 dex). How¬ 
ever, a wide range of metallicities throu ghou t the stream 
are identified by both Ghou et al. ( 2007| ) and Casey et al. 


(2012a). This broad range ot meta l licitie s supports the 
theory of Martmez-Delgado et al. (2004), who suggest 
that the progenitor to the Sgr dwart and stream shed 
its layers (over which there must have been an intrin¬ 
sic MDF gradient) successively in its orbit around the 
Milky Way. This could indicate a mixed orbit al phase in 
the leading arm as described by the model ofjLaw et al.| 

([M^. 

i’he core region of Sgr is predicted to be more metal 
rich than stars in the stream as the stream stars would 
be represented b y debris lost some 3.5 orbits 2.5 — 3 
Gyr) ago (e .g.. Law et al. 2005 2009 Law & Ma¬ 


jewski 2010). Though the mean metallicity over the 
leading arm has been found t o decr ease with distance 
from the core by Chou et al. (2007), the trailing arm 
was not mapped, and the wide range of metallicities 
present in the data prevented a strong detection of a 
gradient. [Keller et al.j (2010) show a significant gra¬ 
dient in the stream at —(2.4 ± 0.3) x 10“^ dex/degree 
(—(9.4 ± 1.1) X 10“^ dex/k pc) from the main b ody of 
Sgr. This is interpreted by Keller et al. (2010) as in¬ 
dicative of a similar gradient in the progenitor of the Sgr 
dwarf and stream system. 


1.1. Target Selection 
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Fig. 1. — Sgr dwarf and stream regions observed are shown as green circles; each observation is 2 degrees in diameter. The Sgr dwarf core 
is located in the center yellow circle and the central rectangular block makes up the ‘central r egion’ FOl. The magenta circles are the first 
obser vations taken, us ed to show no rotation in the core of Sagittarius (|Peharrubia et al.|2011|l. Background image constructed from online 
data (|Mellinger|2009ll and uses a Mercator type projection. We show tlie coordinate grid tor Galactic 1, b, positions where zero marks the 
center ot tne Milky Way. 


In this paper, we present a new selection of probable 
Sgr member stars in the stream using a statistical se¬ 
lection method, and show their metallicity distribution. 
Individual star velocities of the stream will be discussed 
in future stream analysis paper (R. Ibata et al. in prep). 
Our mapping of the Sgr stream includes both leading 
and trailing arm stars as shown in Figure We examine 
phase mixed areas of the stream as well as gradients in 
metallicity along the stream (away from the core). We in¬ 
clude information for known features (overdensities and 
streams) that are near the Sgr stream at va rious points 


on th e sky, namely the Sextans Dwarf (e.g.,|lTwin et al. 


19901, the Virgo Overdensity (e.g., Kellerl 


I20W 


and the 


globular cluster Palomar 5 (e.g., Odenkirchen et ai. 20021 
for comparison. 

A brief overview of the data reduction and analysis is 
outlined in Sections and and the selection of stream 
candidates is given in Section The core sample (the 
center region shown in the outlined box of Figure R is 
hereafter described as ‘FOl’, the remaining stream fields 
shown are divided into 29 regions forming F02 to F30 as 
shown later in Section [6T1 

In Section we compare our candidate Sgr star prop¬ 
erties to points rendered by models of the stream. The 
resulting velocity and metallicity distributions are then 
presented in Section and compared to the model out¬ 
put from a triaxial hmo. We dis cuss a possible detection 
of the Sextans dwarf in Section 16.21 and summarize our 
conclusions in Section [3 

Candidate Sgr member stars were selected using two 
different methods for the two main observing programs 
that were carried out. Firstly we selected stars in the 
stream as part of a joint observing campaign, and later 
the main core area was observed as part of a second, 
more detailed campaign to look for evidence of rotation 


in the core. 

For stars in the stream we based our selection on pho¬ 
tometry from the 2MASS catalog, supplemented (where 
applicable) with optical photometry from the Canada- 
France Hawaii Telescope (CFHT). Candidate Sgr stream 
stars with 2MASS photometry were selected to reside 
within the color and magnitude confines of the following 
box: J - K = (0.95, 1.35, l.IO, 0.95), K = (13.0, lO.O, 
10.0, 11.125). Targets which possessed CFHT photom¬ 
etry were selected as potential members of the follow¬ 
ing broad evolutionary phases: blue horizontal branch 
(BHB), main-sequence (MS), red clump (RC), and red 
giant branch (RGB) stars. The resulting pointings con¬ 
tained stars spanning the magnitude range 12 < V < 18 
(corresponding to a g= 20 limit). 

The second main observing campaign that was carried 
out was on the core region of the Sgr dwarf. For this 
campaign we observed stars within about ~ 6° of the 
center of the Sgr dwarf. The targets were selected based 
on the following criteria for extinction corrected J and 
K magnitudes from 2MASS: 9 < AT < 13 and J — K > 
(20 — Ar)/90, this gives a resulting range in brightness 
between 12 < V < 15. This selection c orresponds to the 
region chosen by Majewski et al. (2004), and is known to 
display a clear trend in velocities across orbital longitude. 

2. OBSERVATIONS AND DATA REDUCTION 

Once selected, can didates were subse quently targeted 
with the AAOmega (Sharp et al. 2006) fibre-fed multi¬ 
object spectrograph on the 3.9m Anglo-Australian Tele¬ 
scope (AAT). We utilize the 1500V and 1700D gratings of 
the AAOmega instrument with resolutions of R= 8000 
and 10000 in the blue and red respectively (the wave¬ 
length range is 4250-6000 A in the blue and 8450-9000 
A in the red)(Sharp et al.||2006). This enables us to tar- 
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get the spectral regions around the Mg I 5170 and Ca II 
8500 A triplets. The AAOmega spectrograph possesses 
392 fibers available for targets in a single pointing. All 
pointings had at least 25 fibers assigned to sky positions. 
Target stars were observed as allowed using the config¬ 
ure program to assign fibers to stars. The observations 
are shown in Figure [l] The data for both the stream 
and core observing programs were obtained on several 
observing runs spanning 2009-2011. 

The spectral reduction was undertaken using the 
.2d/dr|^ data reduction program of the Australian Astro¬ 
nomical Observatory (Heald 2007). This performs bias 
correction, flat-fielding and optimal fiber extraction. Arc 
lamp exposures allowed for wavelength calibration be¬ 
tween each set of observations. Typically, a pointing was 
observed in sets of three exposures, each with an expo¬ 
sure time of 1200 seconds. Cosmic ray reduction was 
then achieved by median combining the multiple object 
frames. We measure an elimination of ~90% of the con¬ 
tribution for all sky lines in the 2dfdr reduction. 

We use a modifi ed version of the Radial Velocity Ex¬ 


periment (RAVE) (Steinm etz et al. 2006 Zwitter et al. 
2008| Siebert et al. |201l[) pipeline to determine metal- 


licity, surface gravity (iog(g)), and effective temperature 
(Teff) for our target stars. We estimate the velocities 


from a template fit on the Ca II 8500 A triplet using the 
IRAl|^ task fxcor, which we then converted into helio¬ 
centric values (Vhei) using the IRAF task rvcorrect. 

After the heliocentric velocity conversion we transform 
the spatial coordinates into Galactic I, b and convert the 
heliocentric radial velocity into a Galactocentric radial 
velocity (Vgsr) using: 


Vgsr = Vhei + 9 cos(l) cos(b) -I- 232 sin(l) cos(b) -I- 7sin(b) 

( 1 ), 

which has units of km s ^ as described in iBraun fc Bur-1 

(119^1. 


ton 


3. ANALYSIS 

The modified RAVE pipeline metallicities are total 
metallicities, referred to hereafter as [m/H]. The abun¬ 
dances are derived from analysis of the wavelength region 
around the Ca II 8500 A triplet using a library of spec¬ 
tral templates. Once an input heliocentric correction is 
given for the observation in question the RAVE pipeline 
derives a temperature estimate and then iteratively de¬ 
termines the best spectral template matciH The weights 
of the bes t match are determi ned by a 3c routine de¬ 
scribed in Zwitter et al. (2008). Conservative estimates 
of errors for a spectrum with an average signal-to-noise 
ratio (S/N ~ 40) are 400 K in te mperature, 0.5 dex i n 
gravity, and_0.2 dex in metallicity (Zwitter et al.||2008). 

In Figure!^ we compare [m/H] metallicity values from 
the modifieoRAVE pipeline with the literature values for 
several stellar clusters: 47 Tucanae, NGC 288, M30, M2, 
Melotte 66, and NGC 1904 from the Harris Catalogue 
(2010 edition; Harris 1996), as listed in Table The 


^ http://www.aao.gov.au/AA0/2df/aaomega 
^ The IRAF software package is distributed by the National Op¬ 
tical Astronomy Observatories, which are operated by the Associ¬ 
ation of Universities for Research in Astronomy, Inc., under coop¬ 
erative agreement with the National Science Foundation. 

® A. Siebert, private communication. 



Fig. 2.— The [m/H] calibration using six reduced calibration 
clusters. The solid line shows one to one correspondence. Each box 
corresponds to one of the clusters given in Table The horizontal 
line in each box is the median for the data measure on that cluster, 
the 25th percentile is given by the upper edge of the box, the 75th 
percentile is given by the lower edge of the box and the maximum 
and minimum are at the ends of the whiskers. The plus symbols 
give the N sigma outliers in the distribution. 


clusters were observed with the AAT using the same in¬ 
strument settings as the Sgr stream targets. The observa¬ 
tions of these clusters were reduced and processed in the 
same way as the Sgr data. The results of the comparison 
between the literature metallicities and those obtained 
with the modified RAVE pipeline are shown in Figure]^ 
The difference between the calculated median shown in 
Figure and the literature value of the mean metallicity 
for the cluster is given by = [Fe/H]iit - [m/H]mean, 
as shown in Table T] In all cases the difference Am <0.31 
dex. As four of the six clusters in Table [T] have a mean 
greater than the median, there is a tendency in this data 
towards a positive skew (greater metallicity values). Due 
to this skew there is closer correspondence between the 
median [m/H]med and the literature mean value [Fe/H]iit. 
We treat [m/H] as equivalent to [Fe/H] and will use the 
[m/H] as determined by the modified RAVE pipeline as 
metallicity for the rest of this paper. 


TABLE 1 

The Calibration Clusters. 


Cluster 

Num 

[Fe/H] lit 
dex 

[ni/H]nie(j 

dex 

[m/Hjniean 

dex 

Am 

dex 

47 Tucanae 

67 

-0.72 

-0.82 

-0.83 

0.11 

NGC 288 

39 

-1.32 

-1.29 

-1.27 

-0.05 

M 30 

25 

-2.27 

-2.15 

-2.06 

-0.21 

M 2 

17 

-1.65 

-1.48 

-1.34 

-0.31 

Melotte 66 

17 

-0.51 

-0.55 

-0.55 

-0.04® 

NGC 1904 

22 

-1.60 

-1.57 

-1.54 

-0.06 


Note. — Calibrati on clusters used from Harris Catalogue (2010 
edition) (|Harris|1996|l for globular clusters. The Web version of the 
database tor Galactic Open Clusters known as BDA (WEBDA) is 
used for the open cluster Melotte 66. [Fe/H]iit gives the literature 
values for the cluster and Am is the difference between [Fe/Hjia 
and [m/H]mean. The median is shown in the box plot of FigurelH 
The Num column gives the stars observed for each cluster. 
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The uncertainties for velocity measurements are 
known, from the AAOmega spectrograph used to take 
the observations, to be 3 to 5 km s“^ due to the spectral 
resolution of the 1700D AAOmega grating. To check the 
robustness of the velocities, and obtain a more precise 
error estimate, we compare the velocities of 21 different 
stars which were observed 3-4 times. We adopt the aver¬ 
age error ojn ~ 2.7 km s“^ as the error for the velocity 
measurements, Vgrr = 2.7 km s“^. This average error 
will only be used for the purpose of draw ing the gener¬ 
alized Gaussian histograms in Section [TT| 

4. SAGITTARIUS STREAM SELECTION 



Fig. 3. — The cuts in log(g) and Tgff for FOl. The filled circles 
represent t he selected stars. The Dartmouth isochrones (|Dotter| 
|et al.|2008|l shown are a 10 Gyr isochrone at fFe/H] = 0.0, as the 
dasned magenta line, and a 10 Gyr isochrone at [Fe/H] = —0.5, as 
the solid yellow line. Open circles are the stars before selection. 


Using the modified RAVE pipeline we exclude targets 
with null or anomalous Ca II 8500 A triplet line mea¬ 
surements. The RAVE correlation coefficient Rcoef deter¬ 
mines the goodness of the fit between the star spectrum 
and the corresponding template through a st andard cross 
correlation procedure ( Tonry fc Davis|l979 ). Higher val¬ 
ues of Rcoef indicate a better ht, but the value of the 
cutoff depends on the effective temperature of the stars 
involved. For example, hotter stars will have intrinsically 
lower Rcoef values. For our data we use the quality se¬ 
lection on the RAVE correlation coefficient Rcoef A 15, 
and visually inspect the residual spectrum (the difference 
between the spectrum and the RAVE model spectrum). 

As mentioned previously, the Sgr stream observations 
are divided into 29 regions (locations are given in Figure 
[^. The K and M-giant stars are selected from RAVE 
mg(g) and Tgff values. 

To separate only the K and M-giants from the data 
we select log(g) and Tgff as follows: log(g) < 2.3 and 
3500 <Teff< 5000. From previous studies of RAVE data 


® The WEBDA metallicity listed is from |Friel & Janes| |1993| 
but all Harris catalog metallicities incorporate multiple literature 
values as found in the Harris bibliography tables. 


we adopt the standard RA VE cutoff, excludin g metallic¬ 
ities of [Fe/H] < —1.5 dex |Siebert et al [2011 ), although 
our data calibration is consistent with a linear calibration 
down to [Fe/H] ~ —2.2 in Figure]^ 

The selection of the FOl core region is shown in Fig¬ 
ure [^ This high density region is the best illustration 
of our method and Figure]^ shows the FOl Sgr selection 
compared to Dartmouth 10 Gyr isochron es for [Fe/H] 


= 0 ( dashed) and —0.5 (solid) metallicities (Dotter et al. 
20081. The choice of 10 Gyr isochrones is motivated 


by previous fitting showing an age range o f approxi- 


et al. 120031 Hellazzini et ai.|2006l 

Law & Maiewski|2U10t 

Niederste-Ostholt et al.||2010 G 

asey et al. 2012a|). In 


9,086 candidate (open circles) stars. Applying the same 
cuts in Teff and log(g) for the stream fields returns 809 
candidates, giving a total of 7,036 candidate Sgr stars. 
The pileup of stars seen near Teft ~ 3500 K in Figure [^ 
is an artifact of the modified RAVE pipeline, and they 
are excluded based on their low RAVE correlation coef¬ 
ficients. 

For coordinates in this paper we adopt the Sgr coordi¬ 


nate system A, B as defined by Majewski et al. (|2003 ). 
This system is oriented such that the equator is defined 
by the Sgr stream midplane, with the origin of A cen¬ 
tered on the Sgr core. This flattens the stream along the 
sky, enhancing the visibility of features in the stream. 


4.1. Membership Likelihood 

After excluding all but the most likely candidates, it is 
still possible that the stars we have selected are members 
of our own Milky Way halo rather than being the desired 
Sgr stream stars. To identify Sgr stars in our sample 
we look for velocity substructure that deviates from the 
Galactic halo. 

We use generalized Gaussian histograms to compare 
the data to the smooth Galactic halo and to define a 
selection criteria for likely Sgr members from the popu¬ 
lation of stars in a field. For each field, we consider the 
range of velocities between Vqsr = —400 and -1-400 km 
s“^. Dividing that range into 100 bins we calculate the 
distribution: 

Dfieid = E D[i] = 1 (2) 

y ITTU 


over all bins, where i is each measurement Vgsr for 
each star in the field in question. The value V is 
given by each measurement of Galactocentric velocity 
(Vgsr) in the field and N is the number of stars in the 
field. The value of a is the error associated with the 
velocity measurements. The candidate distribution for 
each field is Dgeid = Dpoi) ■•••Df3o- We use the value of 
Verr = 2.7 km s“^ as defined in Sectionjsjand V = Vgsr- 
This produces the red velocity distributions shown in 
Figure]^ 

To compare our observations to the background of the 
smooth Galactic halo we use the velocity di stribution as 


predi cted by the Besangon Galaxy model (Robin et al. 
20031. Although there are other models available, to 
make our comparison we only need a general location of 
the halo in velocity space. We find the canonical Be¬ 
sangon model is well suited to our purposes. We calcu- 
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Fig. 4.— The likelihood calculation for two regions, FOl and F15. 
The red thick solid line is the generalized Gaussian histogram for 
our data in the respective regions. The scaled counts axis gives 
the amplitude of the normalized distributions where we have set 
the distribution sampling to be identical to the data. The m edium 
green solid line is the Besangon model ||Robin et al.|2003|l. The 
thin light blue colored lin e is the likelihood v a lue and 

the thin dark line is the |Law & Majewski| ||2010|l model in the 
same region. The top green dasned line marks tire confidence level 
(Cbea) cutoff for membership to the Besangon distribution (the 
smooth Galactic halo). 


late the Besangon generalized Gaussian histogram, Db, 
for each region we observe where the value of V = Vvbc, 
and Vvbc is Vqsr predicted by the model. When us¬ 
ing the Besangon Galaxy model we select each region 
based on the center coordinates of our fields FOl to F30. 
We select stellar types to include M and K giants, and 
bright giants with the default parameters to represent the 
smooth Galactic halo. When we apply Fquationj^to our 
Besangon selection this results in Db, the green curve in 
Figure]^ This distribution is centered at the same loca¬ 
tion on the sky as a corresponding Sgr observation. For 
each Dgeid there is a corresponding Db produced. 

To calculate whether the star is potentially a member 
of Sgr we compare Db with the Dgeid- We also include in 
Figure the distribution from the triaxial model for the 
Sgr Stream (|Law et al.||2009| |Law & Majewski||2010|) as 
the black line, 'i’he black distribution is drawn from the 
model points corresponding to the spatial region spanned 
(in degrees) by fields 1-30 respectively. We show FOl 


and FI5 as examples of the output in Figure The 
black line distribution is simply the generalized Gaussian 
histogram of the model points which lie within the field 
in question, i.e., for FOl this would be the model points 
within ~ 6° of the Sgr core. 

The Dfieid and Db distributions are used to calculate 
the likelihood of a given star in a region not being in the 
smooth halo. To establish limits on the likelihood we set 
stars with Dfleid < Db to have a PNotHaio = 0 and like¬ 
wise if Db < 0.001 we let PNotHaio = 1- In the region 
where Dgeid ^ Db then we need to consider the differ¬ 
ence in the distributions. To keep an accurate scale of 
likelihood between zero and one we will set a confidence 
level which can be used as a cutoff in the region where 
Dfieid > Db. 

For each field we have an appropriately large sample 
of the Besangon model for that area. To determine the 
distribution for each Besangon field (one corresponding 
to each of the 30 data fields) we draw the Besangon dis¬ 
tribution 1,000 times using the same sample number of 
points as are present in the data (the red distribution). 

Setting the confidence level at G=0.996 corresponds 
to the level where 99.6% of the data in the sample 
lies in 1,000 iterations of sampling the Besangon model. 
The location of this confidence level gives us the cut¬ 
off Cbes = G(0.996). We choose this confidence level as 
it roughly corresponds to a 3cr cutoff for normally dis¬ 
tributed data (but we do not assume anything about the 
distribution a priori). The Cbes cutoff is drawn as the 
dotted green line in Figure]^ 

All data in the red distribution which exceed Cbes 
(from the Besangon model) are defined to be likely 
Sgr stars (rather than part of the smooth halo), so if 
Dfieid > Cbes then PNotHaio = 1- If Dfieid > Db and 
Dfieid ^ Cbes USe. 


PNotHaio 


Dfieid Db 

Db _ Cbes 


(3) 


to calculate the likelihood. Equation statistically or¬ 
ders each data point between the limits of zero to one 
using the maximum and minimum defined by the conh- 
dence level Cbes = C(0.996). This gives us a smooth lin¬ 
ear trend between the PnotHaio = 0 (PO) and PNotHaio = 
1 limits. Hereafter, all cases of PNotHaio = 1 are collec¬ 
tively referred to as PI events. The total PNotHaio is 
shown (the light blue curve) with the scaled normalized 
data from the Besangon model, Db, and the Sgr data, 
Dfieid, as shown in Figure]^ For each region Dgeid we 
then have some set of stars with velocity V that are 
included in the PNotHaio = 1 subset and qualify as PI 
events. Combining results from all of our regions we find 
a total of 5513 core FOl candidates and 106 stream can¬ 
didates, a total of 5619 possible Sgr members. Although 
we assume Sgr membership this may not be th e case for 
all objects, as discussed further in Section 


5. COMPARISON OF MODELS TO THE SAGITTARIUS 
STREAM 

Here we compare four models to the PI selection of 
the Sgr stream. These four models ^ are generated with 
triaxial, spherical, prolate, and oblate Milky Way dark 


^ The data for these models are obtained from the website of 
David Law at: http://www.astro.virgiiiia.edu/ srm4n/Sgr/ 
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Fig. 5.— Comparison between our PI data (di amonds) and the prolate, spherical and oblate models from|Law et al.| (|2005||, as well 
as the triaxial model from|Law & Majewski| (|2010||. Both leading and trailing arm PI debris kinematically overlap with tne models. The 
prolate and triaxial models provides a better correspondence to the cluster of PI points at A 2 :: 250°. 


matter haloes, respectively. T he spherical, pro late and 
oblate models were discussed in Law et ah ( 2005 ) w hereas 
the more recent triaxial mod el is trom Eawetah (20091 
and Law & Majewski (2010). 

Mapping the T1 stars returned by our selection tech¬ 
nique, we find that the leading arm of Sgr (known to be 
particularly sensitive both kinematically and spa tially to 
the shape of th e Milky Way dark matter halo; [Law fc 


Majewskij2010 1 overlaps with the predictions of both the 


triaxial and prolate models in velocity space, as shown 
in Figure All four models show agreement in veloc¬ 
ity as they pass through or near PI points but only 
the prolate and triaxial models show agreement in the 
A, Vgsr — (250°, —200 km s“^) region where we have a 
large number of high likelihood PI Sgr stars. 

Though it is known that the triaxial model does not re¬ 
produce the bifurcation in the Sgr stream we consider it 
to be the be st of these four optio ns. A potential solution 
proposed by Casey et al. (2012a) is that this bifurcation 
may have resulted trom a kinernatic disruption of some 
kind, but if the bifurcation was caused by some more re¬ 
cent event then there would be no need to account for the 
bifurcation feature in formational models of the stream. 
However, there is still no model, disruptive or otherwise. 


that can account for all the observed Sgr stream proper¬ 
ties. 

The triaxial model seems to match both leading and 
trailing arm observations kinematically, and as it is the 
most recent of the models we are considering, we use 
it as a point of comparison for our data in the stream. 
We show the full set of PI coordinates and velocities as 
well as other observations from the literature with the 
predictions of the triaxial model in Figure]^ 

6. DISCUSSION 


As shown in Figure we have a range of calculated 
likelihood values that go from zero to one. The highest 
likelihood objects are PNotHaio = 1 (PI) stars, which give 
us a total of 5,619 very high likelihood Sgr members, 106 
of which are in the stream (i.e., not in FOl). The PI 
selection includes only those stars which lay on or outside 
of the Cbes dotted blue line shown in Figure® 

We find that 100% of the PI stars selectedin the FOl 
region lie within the known velocity range for the Sgr 
dwarf, i.e., within a range of heliocentric velocities be¬ 
tween 100 and 200 km s“^. This range is chosen based 
on th e radial velocity prof ile of the core to within 3cr 
(from ^llazzini et al.|200'8 Vcore = 139.4± 0.6 km s“^). 


Including all Sgr candidates for FOl we find 93.3% of the 
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Fig. 6.— The distribution of the highest likelihood Sgr stars. PI stars are shown as large red diamonds, the Virgo Overdensity locations 
are shown as black pentagons, Pal5 is shown as a thin blue diamond. Top is the field nu mbering used in the likel ihood calculations, starting 
at the core (FOl) a nd going t o stream field 30. Fo r comparison we show the data fromi Maiewski et al.lll20 04ll as squares, the data from 
iKeller et al.| ||201 0|| as circles, iMonaco et al.11[200^ as upwards triangles, and data from luhou et akjll ^UuYl) as downwards triangles. The 
tnaxial mogel ot |L;aw &: Majewski] J^LUl is given by the very small dots (online data from models or |Paw et al.||2009[ [Law &: Majews^ 
|2010| . The trailing arm is shown in light blue and starts at a velocity of about 180 km s“^ at A = 0 and goes down to —150 km s“^ at 
360 degrees, whereas the leading arm is shown in yellow and starts at —150 km s“^ at A = 0 and goes up to 180 km s“^. 


stars are within the velocity range for Sgr. The PI se¬ 
lection then gives us zero halo objects in the well known 
Sgr core region, without assuming anything beyond that 
they are not part of the smooth Milky Way halo. We use 
these PI stars to trace likely Sgr members in the stream 
where the velocity distribution is less well known. The 
spatial and velocity distribut ion of Sgr PI dat a is co m- 
pared to the triaxial model of Law & Majewski (2010) in 
Figure]^ We additionally show the positions and veloci¬ 


ties observed for the stream from Majew ski et al.|(|2004 
Chou et ^ (2007), Monaco et al. (2007), and |Keller et ah 


(20101 tor comparison. _ 

'There are 106 PI stream members shown in Figure 
The PI data indicate that the distribution of Sgr stars in 
the core and stream approximat ely agree with the predic¬ 
tions from the triaxial model of Law & Majewski (2010). 
In both leading and trailing arm regions of the model 
there is a general consistency in velocity space. The dis- 
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crepancy between data and models could be due to prob¬ 
lems in the models or with the small number statistics, 
i.e., the low density of PI stars in the stream. A list of 
the selected PI stars is given in Table 

6.1. The Stream Metallicity Distribution 
We note agreement in velocit y space for the PI stream 


from d i fferent evolutionary phase s is given in Law fc Ma-| 
lewski ( 2010 ). As mentioned by Chou et al. (2007p and 
Keller et al. (2010|, the longer trailing arm yields better 


stars and the triaxial model ( Law fc Majewski |2010p , 
listed as ‘Law’ in the bottom panel ol h'lgurej^ 1 hough 
there are a few interesting discrepancies (these may be 
due to other unknown overdensities or kinematic sub¬ 
structure, as discussed in Section 
that the majority of PI stars wil 


6.2), the expectation is 
be Sgr members. 

The largest discrepancies between the data and the 
triaxial model lie in the region A ~ 150° — 250° shown in 
Figure and Figure The triaxial model trailing arm 
starts at a velocity of about 180 km s“^ at A = 0 and 
goes down to a velocity of about —150 km s“^ at 360 
degrees whereas the leading arm starts at —150 km s“^ 
at A = 0 and goes up to 180 km s“^. The areas with the 
highest density of PI stars show a good correspondence 
with the model but the mis-matches indicate that there is 
more work to be done in particular in fitting the trailing 
arm. Towards that end we investigate the metallicity 
distribution from the core and across the stream. 

We find the average metallicity of core (FOl) PI stars 
to be —0.59 dex with a dispersion of ct ~ 0.34 dex. This is 


sured mean 

of [Fe/Hl~ —0.5 dex (Monaco et al. 

2005 

Keller et al. 

2010| Cacciari et al. 

2002), and we expect 


indeed may even go into the extremely metal poor range, 
i.e. having metallicities of [Fe/H] < —3.0. 

Without selecting purely for M-giants we have about 
30 stars with [m/H] < —2.0 dex in our pre-processed FOl 
data. These stars have velocities and positions consistent 
with the central region of the Sgr dwarf (within ~ 6° of 
the core), and these intrig uing objects are p art of an on¬ 
going follow-up campaign ( Hyde et al.|2012| ). The RAVE 


errors on many of these measurements do not pass the 
selection criteria for this paper but the spectra are sorted 
via an individual splot IRAF measurement of the Ca 11 
8500 A triplet spectral lines. 

If our 30 potentially metal poor Sgr objects can be con¬ 
firmed with high resolution spectroscopy, this number of 
metal poor objects i n FOl would support the findings of 


Casey et al. (2012a), who indicate that the distribution 
m the core of Sgr may be more metal poor than was 
previously thought. As stated previously, for this publi¬ 
cation we disregard any poor fits to the RAVE templates. 
This means many of the very lowest metallicity stars fall 
out of the reduction, as they are not measurable by the 
automated routines we used. 

To investigate the leading and trailing arm of the 
stream we select stream stars that are coincident with 
the triaxial model in the A, B, and Vgsr coordinates. 
We also require that the stars not be in an overlapping 
region, i.e., areas of the model where both leading and 
trailing stars are predicted to lie. In EigurelTl we show the 
bins containing only PI points. Law et ahj 
that along the leading arm of t. 


Law et akf ( |2005| ) note 
le stream, stars which 
are lost from Sgr in different orbits around the Milky 
Way can overlap in orbital phase position. Possible ef¬ 
fects of angular phase smearing and populations of stars 


energy sorting of the debris and can be more cleanly iso¬ 
lated from background or mixed populations. To select 
a clean set for leading and trailing arm populations we 
choose only stars which do not occupy a mixed phase 
position and create the leading and trailing population 
as shown in Eigurej^ 

In Figure we create 2 bins for the leading arm 
(green) and 3 bins for the trailing arm (light blue) 
along 360 degrees of A. These bins are defined to 
maximize the number of stars included as well as the 
spatial resolution as shown in Figure The five 

stream bins are located at A ~ [301, 246,142,106,14] 
with a range of Amin — Amax — [30,17,17,18,1] and 
[13,10,11,24,16] stars in each bin respectively. The dis¬ 
tributions in metallicity in these bins are shown in Fig¬ 
ure [^ with the median of the bin given by a horizontal 
line. The mean metallicity of the five stream bins are 
[-0.27, -0.60, -0.97, -0.63, -0.64] dex, with MDF dis¬ 
persions of [0.55,0.52,0.26,0.24,0.33] dex corresponding 
to the A coordinates above. The core FOl region has a 
mean metallicity of —0.59 dex and a dispersion of 0.34 
dex from our sample of 5513 PI core stars. The PI ob¬ 
jects which do not exclusively lie within a leading or trail¬ 
ing arm portion of the model are not included in the bins. 

We find that the mean metallicity starting from —0.59 
dex in the core, drops from —0.63 dex to —0.97 dex as 
A increases along the trailing arm in Figure (in the 
positive direction). In the leading arm, however, there is 
little or no trend with A for the metallicity. The higher 
variation in metallicity values in the leading arm may 
be preventin g any discernment of trends. Mapping of 
the arms by Chou et al. (2007) indicates that leading 
arm stars seem to be more metal poor than t he core. In 
our data as well as that of Chou et al. (2007|) the differ¬ 
ence in the mean metallicities of the sarnples is less than 
their MDF dispersions (0.31 and 0.33 dex respectively), 
indicating that more detections will be needed before a 
gradient can be considered to be well established. 

Kinemati c s of red horizontal branch stars from 
Shi et al. ( [2012 ) suggest a metallicity gradient of 
— (1.8 ±0.3) X 10“^ dex deg“^ in the trailing arm and 
a smaller gradient of —(1.5 ± 0.4) x 10“^ dex deg“^ in 
the leading arm. Although these methods are different 
from the ones used in this paper, such a gradient does 
agree with theories for the formation of the dwarf galaxy 
where a gradient is initially present. In particular, the 
more metal poor stars sh own in ou r traili ng arm data 
agree with the scenario of Shi et al. (2012) wherein the 
Sgr dwarf galaxy would have been formed with a more 
metal rich core and would have been surrounded by older 
and more metal poor stars. Those more metal poor stars 
would then have been stripped first and found further 
al ong the stream, c reating the gradient we see. 


Chou et al. (2007) find that the MDF of the Sgr stream 
changes frorn a median of [Fe/H]~ —0.4 dex to about 
— 1.1 dex over a leading arm length, but our values do 
not drop off at that level, staying in the range of [m/H]~ 
—0.61 dex to about —0.27 dex. Though the metallicity 
distribution in leading arm stars is known to be quite 
broad we find that our leading arm stars seem to lack a 
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Fig. 7. — MDF evolution with A for the leading (two negative A boxes) and trailing (three positive A boxes) Sgr arms. The Sgr core 
(FOl) region is shown by the single red boxed area at A = 0. The black line in each box is the median value of metallicity for the area, 
the 25th percentile of the data is given by the upper edge of the box, the 75th percentile is given by the lower edge of the box and the 
maximum a nd minimum are a t the ends of the whi skers. Our Sgr point s included are sho wn as yellow diam onds, comparison data are 
shown from [filler et aT]||2010|| as small blue circles, [Monaco et al.|([2007^ as diamonds, andjChou et al.|||2007[l as pentagons. 


decrease to metal poor, and do not seem to show more 
metal poor s tars t han the trailing arm as mentioned by 


Chou et al. (2007). The broad range of metallici ties in 
the leading a rm has additionally been noted by |Casey| 


et al. (2012a). The median metallicity of the Galactic 
thick disk is ^ —0.7 dex, and, as noted by |Chou et al.| 
(2007), should not have a large impact on MDF trends 
across the stream. The trend to lowe r metallicities in the 
leading arm of the stream found by Chou et al. (2007) 
may be related to their slightly more metal rich value tor 
the Sgr core. 

The comparison of the Sgr core region with the lead¬ 
ing arm does not show the decline in median metallicity 
which would be repr esented by debri s lost some 3.5 orbits 
2.5 — 3 Gyr) ago (Law et al. 2005), but we do find sub¬ 
stantial variation in the iVlDl' along the strea m. The ob¬ 
served MDF variatio n supports the theory of [Martinez- 1 


Delgado et al. (2004), who suggest the satellite shed suc¬ 
cessive layers in its orbit, over which there must have 
been an intrinsic MDF gradient. The disagreement in 
the direction of the trend between our data and that of 


Chou et al.|(2007l indicates that perhaps the variation 
may be wider, and the gradient less, than either of our 
samples suggest. This scenario would support models in 
which a rapid change in the binding energy of Sgr oc¬ 
curred over the past several gigayears, providing a large 
net metallicity variation, but a shallow gradient. The 
sudden change of state could have been caused by some 
dramatic event, which may als o provide the source of 
kinematic disruption invoked by Casey et al. (2012a) to 
account for the bifurcation seen in the streani. 

We additionally find a large MDF variation in the trail¬ 
ing arm sample, as mentioned above. If this region is 
truly less susceptible to overlap in orbital phase position 
then this would indicate evolution in the stream from a 


mean of -0.63 (near what we find for the core) to -0.97 dex 
(our most distant trailing arm box). This variation is ad¬ 
ditionally on the order of the MDF dispersion, but with 
overall lower dispersion than what is found in the leading 
arm. We then can conclude that although the dispersion 
in the leading arm may be due to mixed orbital phase, 
the trends we find in the leading and trailing arms sup¬ 
port a change in the Sgr MDF with position along the 
stream, and therefore an intrinsic MDF gradient in the 
progenitor to the Sgr dwarf and stream system we see 
today. 


6.2. Overdensities in the Stream 

In Figure]^ we include information for known features 
(overdensities and streams) that are near the Sgr stream 
at various points on the sky. While we don’t target them 
in our observations, there are several which have note¬ 
worthy overlaps. 

Firstly we consider the Vi rgo Overden s ities (VOD). 
Discovered indep e ndently by [V ivas et al.j (|20()T1); |New -| 
berg et al.[([20021);[Keller ([201Up and|Vivas et al.[ (|2U12|), 
the VOD includes three halo substructures: the two 
structures at an Right Ascension of 160 and 180 degrees 
respectively (at distances of 17 and 19 kpc, and with 
radii of 1.3 and 1.5 kpc) and an extended feature at 28 
kpc that covers at least 162 deg^ (the Virgo Equatorial 
Stre am). Perhaps covering 1000 sq deg or more on the 
sky ([Juric et al.|20 08[[Bonaca et al.|2012 ). We use the co- 
ordinates found in Keller } 2010) to plot the spatial points 
and we adopt Vosr = 130±10 km s“^ for the VOD from 


Newberg et al. (2007). This gives us the black pentagons 
shown in k igure^ 

The Palomar 5 (Pal5) globular cluster is located at 
coordinates 1, b ^ (0.85°,45.86°) shown by the light 
blue diamond. The metallicity is [Fe/H]=—1.41 dex from 
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TABLE 2 

The approximate position and velocity of PI overdensities. 


Field 

Stars 

Aavg 

Pavg 

R A 

-‘-^■'^avg 

DECavg 

lm/HJ„„g 


VGSR(avg) 

tTv 

catagory 



degrees 


degrees 


dex 

dex 

(km 8“^) 

(km s-^) 


01 

5513 

154.32 

1.65 

19.0 

-30.8 

-0.59 

0.34 

168.76 

12.15 

Core FOl region of Sgr 

04 

23 

103.05 

-0.22 

02.0 

-01.9 

-0.65 

0.2 

-112.61 

10.90 

Potential Sgr Stream clump 

15 

17 

231.02 

27.52 

10.2 

-01.4 

-1.10,±0.48 

0.65 

79.45 

11.49 

Potential Sextans detection 

23 

3 

300.31 

0.23 

15.0 

-06.9 

0.38 

0.02 

50.71 

5.39 

Potential Sgr Stream clump 


Note. — These groups correspond only to collections of three or more stars within ±25 km s ^ in velocity for an area with a 5 degree 
radius on the sky. The distribution of all points is given in Figure]^ 


the Harris Catalogue (2010 edition; Harris 19961 which 
is only slightly richer than the Sgr core region. Pal5 
has an extremely low velocity dispersion, with a helio¬ 
centric velocity Vhei = —58.7 ± 0.2 km s“^ and a to- 
tal line-of-sight velocity d ispersion of 1.1 ± 0.2 km s“^ 
(Odenkirchen et al. 2002). From Equation this trans¬ 
lates into Vgsr ~ —106.47 km s“^. 

Finally we consider the Sextans dwarf, centered on 
the light blue circle. At coordinates of (RA, Dec) ~ 
(10/il3', -01°36') or (l,b) ~ (243°.5, 42°.3) the Sex- 


tans dwarf is known to have a Galactocentric veloc- 

itv of Vgsr —73 km s ^ (Irwin et al.||l990 

Irwin & 

Hatzidimitriou| 1995 

Karachentsev et al.||2l)U4 

. Sextans 

members have been 

measured to have a metallicity ot 


[Fe/H]= —1.7 ± 0.25 dex (Da Costa et al. 1991) which 
is slightly more metal poor than what we expect for Sgr 
stars. With a diameter of 30 arcminutes on the sky, this 
large and somewhat diffuse structure is similar to the 
type of detection we are looking for in the Sgr stream. 

The structure of Sgr overlaps with these overdensities 
in several regions and it appears to remain well separated 
when metallicity, velocity, and spatial coordinates are 
taken into account, except in the case of the Sextans 
dwarf. The presence of the comparison features above 
led us to search for additional groupings of stars in our 
data. 

We searched the PI data and found 4 regions which 
have stars that are coincident in velocity and spatial co¬ 
ordinates; i.e. within ±25 km s“^ in velocity for an area 
with a 5 degree radius on the sky. These sets or ‘over 
dense’ regions in the data are located in FOl, F04, F15 
and F23 as shown in Table While F01,F04, and F23 
velocities correspond with what is expected from predic¬ 
tions of the triaxial model, the group of ~ 17 stars from 
F15 has a large offset with respect to the triaxial model. 
This detection is easily observed in Figure 13 where com¬ 
parison shows little correspondence with the Milky Way 
smooth halo or the triaxial model for Sagittarius. The 
peaks for the MDF and the velocity distribution of F15 
are given in Figure 

The group of 17 stars in F15 has an average velocity of 
Vgsr = 79.45 km s“^ (where a = 11.49 km s“^) and an 
average metallicity of [m/H]ai;g = —0.83 (where a = 0.65 
dex) for a position centered on A, B ~ (231.02°, 27.52°). 
From Figure the MDF seems to be bimodal; however, 
the three supersolar metallicity values do not reflect a 
grouping in velocity. The coordinates in ([m/H],VGSR) 
are (±0.39,83.79), (±0.5,65.31), and (±0.5,103.94) and 
removing them from our distribution will shift the mean 
metallicity to [m/H]~ —1.10 dex but will have almost no 
effect on the average velocity. 

The FI5 group is near the spatial coordinates of three 
of our overdensities, namely the VOD and the Sex¬ 



Fig. 8.— The top panel plots the MDF for the F15 overden- 
sity, showing a super-solar population as well as a lower metallicity 
group at [m/H]~ —1.10 dex. The bottom panel shows the velocity 
distribution for the F15 overdensity; the higher velocity objects do 
not correspond to the super solar metallicity group. The detection 
of PI stars for F15 in Figure translates to the above ranges in 
[m/H] and Vqsr- 

tans dwarf. However, the velocity measured corresponds 
closely only with Sextans. As mentioned previously, we 
expect a value of approximately 73 km s“^ from the liter¬ 
ature for Sextans, and we measure Vgsr = 79.45 km s“^ 
for F15. This is dissimilar to the other overdensities in 
the region, i.e., the adopted VOD velocity (Vgsr — 130 
km s“^) and the expected velocity of the Sgr stream from 
the triaxial model (which shows peaks at Vgsr — 0 and 
— 150 km s“^ as shown in Figure]^. This potential de¬ 
tection of Sextans in our data provides an additional 
validation for our statistical method and its ability to 
distinguish stellar structures from the smooth Galactic 
halo. The diffuse nature of Sextans provides us with an 
additional object which has an extended nature and suf¬ 
fers from foreground star contamination. That we can 
detect this object with a strong signal lends confidence 
to our other detections as truly being ‘nof part of the 
Milky Way halo population. 

7. CONCLUSIONS 

We introduce a selection technique to separate Sgr 
member stars from the background. As part of this 
method we first select only the best K and M-giant mea¬ 
surements from 24,110 spectroscopic observations. We 
then assign a likelihood of a given star not being part of 
the smooth Galactic halo. This is done by drawing a dis- 
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tribution using Equation and calculating a likelihood 
per star per field using Equationj^ We identify 106 likely 
members of the Sgr stream and nnd three main results: 


1. The Sgr stream Vqsr distribution with A shows 
reasonable agreement with a triaxial model for the 
Milky Ways dark matter halo. We compare the 
kinematics of the observed Sgr stream stars with 
those of extant simulations of the tidal disruption 
of Sgr. From compariso ns with 4 models we id en- 
tify the triaxial model of Law & Majewski (20101 as 
a reasonable fit to the data and find general agree- 
ment across the Sgr stream. 


2. The Sgr stream MDF variation with A has yielded 
several interesting results. We find that the trail¬ 
ing arm mean metallicity seems to become more 
negative further away from the Sgr core. We mea¬ 
sure a decrease in mean metallicity from —0.59 dex 
in the core to —0.97 dex with increasing A. This 
metallicity gradient, in which material further from 
the Sagittarius core is more metal-poor, is consis¬ 
tent with the scenario of tidal disruption from a 
progenitor dwarf galaxy that possessed an internal 
metallicity gradient. In the leading arm, however, 
we have larger variations in metallicity values, lack¬ 
ing a clear trend; this larger range of values may 


be due to a mixing of orbital phases in the leading 
arm. 


3. A search for overdensities finds three groups con¬ 
sistent with the triaxial model (FOl, F04, F23) 
and one potential detection of the Sextans dwarf 
at F15. We report on this new detection, the F15 
overdensity, located at RA, DEC ~ (10.2, —01.4), 
which has peaks in the metallicity distribution at 
[Fe/H] = —1.10 and -1-0.46 dex and an average ve¬ 
locity of Vqsr = 78.48 km s“^. The F15 overden¬ 
sity does not appear to coincide with e ither the tri¬ 
axial model of Law & Majewski (2010), the VOD or 
the smooth Galactic halo as represented by the Be- 
sangon model; but it does coincide in both spatial 
and velocity coordinates with the Sextans dwarf. 


The parameter space we use for this study includes the 
modified RAVE pipeline log(g), Tgff, [m/H] metallicity, 
and velocity. We obtained velocity information through 


fitting the Ca II 8500 A triplet with IRAF software and 
we use the 2MASS photometric colors for initial object 
selection, but for many of our obje cts this parameter 
space could be expanded. The ppmxl { Roeser et al.|2010| 


catalo g, the Sloan Digital Sky Survey (SDSS) ( York et al. 


20001 observat ions, and the ongoin g Skymapper photo- 


metric survey (Casey et al. 2012bI may all be incorpo¬ 
rated in the future. 


TABLE 3 The selected 106 likely Sagittarius stream members with modified 
RAVE pipeline metallicity, surface gravity and effective temperature values. 


Num 

i^Aj2000 

DECj2000 

A 

(degrees) 

B 

(degrees) 

log(g) 

Teff 

(K) 

[m/Hj 

(dex) 

Rcoefi 

1 

01:58:05.38 

-01:45:11.2 

102.62 

-0.58 

1.31 

4249.0 

-0.92 

21.21 

2 

01:58:37.34 

-01:44:19.0 

102.75 

-0.52 

1.18 

4182.0 

-0.49 

28.49 

3 

01:58:37.34 

-01:44:19.0 

102.75 

-0.52 

1.09 

4209.0 

-0.65 

31.0 

4 

01:58:45.53 

-01:43:24.2 

102.78 

-0.52 

1.48 

4220.0 

-0.4 

23.83 

5 

01:58:45.53 

-01:43:24.2 

102.78 

-0.52 

1.34 

4148.0 

-0.62 

22.76 

6 

01:59:10.97 

-02:19:21.7 

102.57 

0.05 

1.66 

4571.0 

-0.55 

36.68 

7 

01:59:50.95 

-02:08:49.6 

102.8 

-0.02 

1.08 

4098.0 

-0.12 

20.4 

8 

02:00:07.30 

-02:23:16.8 

102.74 

0.23 

1.36 

4206.0 

-0.44 

31.64 

9 

02:00:10.75 

-01:39:51.5 

103.12 

-0.39 

1.29 

4470.0 

-0.71 

30.25 

10 

02:00:10.75 

-01:39:51.5 

103.12 

-0.39 

1.66 

4511.0 

-0.9 

36.04 

11 

02:00:33.24 

-01:36:46.4 

103.23 

-0.39 

1.75 

4368.0 

-0.81 

21.76 

12 

02:00:33.24 

-01:36:46.4 

103.23 

-0.39 

1.95 

4330.0 

-0.89 

23.38 

13 

02:00:36.79 

-02:05:40.9 

102.99 

0.03 

1.16 

4378.0 

-0.65 

39.38 

14 

02:00:36.79 

-02:05:40.9 

102.99 

0.03 

0.91 

4352.0 

-0.91 

29.51 

15 

02:01:00.19 

-01:48:00.0 

103.23 

-0.17 

2.2 

4674.0 

-0.67 

24.61 

16 

02:01:00.19 

-01:48:00.0 

103.23 

-0.17 

1.94 

4891.0 

-0.82 

25.6 

17 

02:01:03.12 

-02:02:04.2 

103.12 

0.04 

0.27 

3926.0 

-0.95 

41.44 

18 

02:01:03.12 

-02:02:04.2 

103.12 

0.04 

1.67 

4304.0 

-0.55 

34.22 

19 

02:01:16.68 

-01:34:20.6 

103.4 

-0.33 

1.78 

4279.0 

-0.47 

21.67 

20 

02:02:33.29 

-01:38:29.4 

103.64 

-0.11 

2.22 

4505.0 

-0.49 

29.09 

21 

03:00:12.67 

+05:50:44.5 

119.85 

0.69 

1.62 

3836.0 

-0.65 

39.23 

22 

03:00:20.59 

+06:02:06.7 

119.97 

0.54 

2.15 

4404.0 

-0.13 

15.83 

23 

03:00:40.32 

+06:34:28.6 

120.31 

0.11 

1.84 

4679.0 

-0.25 

27.24 

24 

03:02:50.71 

+06:53:02.0 

120.93 

0.11 

0.82 

4116.0 

-0.99 

34.77 

25 

03:58:20.57 

+13:42:56.9 

136.18 

0.64 

1.11 

4559.0 

-1.5 

39.35 

26 

03:58:45.77 

+13:58:35.4 

136.39 

0.46 

0.54 

3791.0 

-1.17 

43.81 

27 

03:59:42.72 

+14:07:32.5 

136.66 

0.43 

1.33 

4020.0 

-0.59 

19.55 

28 

04:00:24.00 

+14:40:31.1 

137.06 

0.01 

1.51 

4482.0 

-0.77 

37.33 

29 

04:00:58.70 

+13:52:26.4 

136.82 

0.79 

0.68 

3910.0 

-1.02 

36.89 

30 

04:00:58.73 

+13:52:26.8 

136.82 

0.79 

0.63 

3932.0 

-0.93 

35.94 

31 

04:03:19.70 

+14:56:39.1 

137.81 

0.09 

0.26 

3934.0 

-1.08 

27.64 

32 

04:59:05.93 

+21:28:01.6 

152.56 

-0.28 

1.25 

4263.0 

-0.62 

18.15 

33 

04:59:06.00 

+20:30:37.1 

152.2 

0.6 

1.52 

4426.0 

-0.74 

23.61 

34 

05:00:29.88 

+21:20:42.4 

152.82 

-0.05 

0.35 

3923.0 

-1.22 

35.07 

35 

05:00:29.88 

+21:20:42.4 

152.82 

-0.05 

0.73 

3977.0 

-1.01 

33.35 

36 

05:02:03.98 

+21:04:52.0 

153.05 

0.34 

1.29 

4547.0 

0.5 

20.47 

37 

06:59:45.96 

+21:46:16.7 

179.6 

7.27 

2.28 

4353.0 

-0.52 

25.18 

38 

07:00:27.79 

+21:46:32.5 

179.76 

7.29 

0.81 

4421.0 

-1.23 

28.01 


Continued on next page 
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Num 

RAj2000 

DECj2000 

A 

B 

log(g) 

Tefif 

[m/H| 

Rcoeff 

39 

07;00;27.79 

+21:46:32.5 

179.76 

7.29 

1.86 

4836.0 

-0.86 

45.42 

40 

07:03:33.91 

+21:34:03.7 

180.44 

7.61 

1.86 

4797.0 

-1.28 

49.39 

41 

07:59:42.55 

+32:51:33.5 

193.94 

-2.43 

1.37 

4702.0 

-1.65 

24.28 

42 

10:00:04.22 

+18:18:12.2 

222.36 

9.25 

2.11 

4527.0 

0.28 

36.88 

43 

10:00:04.22 

+18:18:12.2 

222.36 

9.25 

2.25 

4571.0 

0.24 

34.09 

44 

10:11:48.29 

-01:37:32.5 

230.65 

27.81 

2.1 

3894.0 

0.5 

16.61 

45 

10:11:52.85 

-01:44:28.3 

230.71 

27.92 

0.84 

3942.0 

0.39 

22.06 

46 

10:11:52.85 

-01:44:28.3 

230.71 

27.92 

2.17 

3997.0 

0.5 

21.46 

47 

10:12:26.21 

-01:38:30.1 

230.83 

27.78 

1.86 

4117.0 

-1.07 

17.5 

48 

10:12:41.59 

-01:32:51.7 

230.87 

27.68 

1.14 

3974.0 

-1.72 

23.26 

49 

10:12:41.59 

-01:32:51.7 

230.87 

27.68 

1.98 

4262.0 

-0.89 

23.36 

50 

10:12:41.83 

-01:45:27.4 

230.93 

27.88 

1.87 

4720.0 

-1.11 

37.92 

51 

10:13:20.74 

-01:31:59.2 

231.04 

27.62 

2.09 

3509.0 

-0.61 

19.4 

52 

10:13:22.94 

-01:22:27.1 

231.0 

27.46 

1.64 

4035.0 

-1.18 

20.94 

53 

10:13:34.68 

-01:01:59.2 

230.94 

27.12 

1.94 

4175.0 

-0.65 

21.71 

54 

10:13:35.30 

-01:17:54.6 

231.03 

27.37 

1.43 

4497.0 

-1.11 

24.19 

55 

10:13:55.10 

-01:17:03.8 

231.11 

27.34 

1.04 

4250.0 

-1.34 

26.38 

56 

10:13:55.10 

-01:17:03.8 

231.11 

27.34 

2.21 

4685.0 

-1.24 

19.92 

57 

10:14:04.30 

-01:13:22.4 

231.14 

27.27 

1.58 

4080.0 

-1.38 

18.51 

58 

10:14:04.68 

-01:23:46.7 

231.19 

27.43 

1.77 

4179.0 

-0.87 

18.15 

59 

10:14:39.86 

-01:10:52.7 

231.28 

27.18 

1.05 

4321.0 

-1.0 

19.66 

60 

10:15:52.08 

-01:33:31.0 

231.73 

27.46 

1.18 

4113.0 

-1.29 

24.13 

61 

10:58:20.90 

+15:34:44.8 

236.66 

7.91 

1.18 

4243.0 

-0.7 

29.94 

62 

10:58:20.90 

+15:34:44.8 

236.66 

7.91 

1.91 

4470.0 

-0.61 

37.86 

63 

11:01:45.26 

+16:00:08.3 

237.29 

7.23 

1.32 

4604.0 

-1.11 

33.88 

64 

11:02:19.32 

+15:31:30.0 

237.58 

7.64 

1.27 

3918.0 

-0.53 

40.89 

65 

11:02:19.34 

+15:31:30.4 

237.58 

7.64 

1.34 

3959.0 

-0.51 

32.07 

66 

11:02:27.36 

+15:35:13.6 

237.59 

7.57 

1.66 

4833.0 

-1.48 

43.44 

67 

12:25:28.20 

+24:20:53.2 

251.97 

-8.25 

1.66 

3677.0 

0.01 

15.43 

68 

12:25:35.88 

+23:18:00.0 

252.48 

-7.33 

2.14 

3559.0 

0.39 

26.38 

69 

12:26:20.54 

+24:31:34.3 

252.07 

-8.5 

1.64 

4756.0 

-0.85 

16.98 

70 

12:27:24.17 

+23:16:01.9 

252.87 

-7.49 

1.63 

4082.0 

-0.48 

19.86 

71 

12:29:41.66 

+24:04:06.6 

252.97 

-8.44 

1.09 

4347.0 

-1.27 

26.03 

72 

12:29:41.66 

+24:04:06.6 

252.97 

-8.44 

0.96 

4311.0 

-1.49 

27.65 

73 

14:12:45.12 

-00:46:50.5 

287.07 

0.91 

2.1 

3893.0 

0.5 

19.36 

74 

14:14:29.23 

-00:31:49.4 

287.31 

0.47 

1.57 

4299.0 

0.03 

25.72 

75 

14:16:22.44 

-01:04:48.0 

288.0 

0.7 

0.42 

3987.0 

-1.38 

39.15 

76 

14:16:34.99 

-00:21:00.4 

287.68 

0.05 

1.98 

3640.0 

0.43 

24.9 

77 

14:16:36.60 

-00:10:30.4 

287.59 

-0.11 

2.03 

4592.0 

-0.01 

20.38 

78 

14:59:10.44 

-06:47:11.4 

300.09 

0.26 

2.17 

3688.0 

0.38 

35.51 

79 

14:59:58.13 

-07:10:29.3 

300.46 

0.5 

2.21 

4801.0 

-0.3 

24.55 

80 

15:18:34.13 

+00:17:22.9 

300.8 

-8.29 

0.95 

3857.0 

-0.72 

33.39 

81 

15:59:00.53 

-14:46:18.8 

316.8 

0.23 

1.37 

4225.0 

-0.62 

22.78 

82 

15:59:31.73 

-14:18:20.5 

316.7 

-0.25 

0.98 

3791.0 

-0.69 

28.23 

83 

16:00:05.66 

-14:39:53.3 

316.99 

0.01 

0.97 

4060.0 

-0.9 

27.88 

84 

16:00:46.25 

-10:08:09.2 

315.08 

-4.1 

1.99 

4315.0 

0.07 

16.32 

85 

16:01:23.09 

-14:28:50.5 

317.18 

-0.3 

1.92 

4214.0 

-0.23 

30.47 

86 

17:00:54.10 

-19:57:31.3 

332.38 

-1.27 

1.56 

4211.0 

-0.7 

32.8 

87 

19:57:50.11 

-33:16:26.0 

13.56 

2.86 

1.32 

4319.0 

-0.66 

23.72 

88 

19:58:16.32 

-32:46:55.6 

13.64 

2.37 

1.8 

4499.0 

-0.2 

25.2 

89 

19:58:17.09 

-33:16:32.5 

13.66 

2.86 

2.07 

4177.0 

-0.33 

15.33 

90 

19:58:39.50 

-32:59:30.1 

13.72 

2.57 

1.84 

4362.0 

-0.67 

17.45 

91 

19:58:50.18 

-33:05:16.1 

13.77 

2.67 

1.39 

4225.0 

-1.23 

16.2 

92 

19:59:02.09 

-33:25:51.6 

13.82 

3.01 

1.95 

4646.0 

-0.98 

52.36 

93 

19:59:48.53 

-33:20:13.6 

13.98 

2.91 

1.67 

4008.0 

-0.31 

23.75 

94 

20:00:08.81 

-32:42:50.8 

14.03 

2.28 

1.8 

4370.0 

-0.68 

33.01 

95 

20:00:34.25 

-32:56:00.2 

14.12 

2.5 

1.95 

4753.0 

-0.63 

29.79 

96 

20:00:48.65 

-33:19:37.2 

14.19 

2.89 

2.25 

4723.0 

-1.03 

30.25 

97 

20:01:06.31 

-32:41:03.8 

14.23 

2.25 

1.89 

4264.0 

-0.6 

27.7 

98 

20:01:08.74 

-32:39:25.6 

14.24 

2.22 

1.61 

4241.0 

-0.45 

20.47 

99 

20:01:12.02 

-32:50:30.8 

14.25 

2.41 

2.16 

4248.0 

0.02 

16.54 

100 

20:01:20.64 

-32:36:18.7 

14.28 

2.17 

2.18 

4363.0 

-0.5 

31.93 

101 

20:01:20.81 

-33:20:58.9 

14.3 

2.91 

1.89 

4586.0 

-0.82 

28.07 

102 

20:02:09.79 

-33:13:50.9 

14.47 

2.79 

1.1 

4368.0 

-1.13 

38.56 

103 

22:00:16.66 

-30:08:30.5 

39.72 

2.29 

0.81 

4300.0 

-1.25 

40.51 

104 

22:00:16.66 

-30:08:30.5 

39.72 

2.29 

1.08 

4299.0 

-1.15 

52.7 

105 

22:00:22.90 

-30:07:50.9 

39.74 

2.28 

1.1 

4513.0 

-1.47 

38.94 

106 

22:00:22.90 

-30:07:50.9 

39.74 

2.28 

1.03 

4295.0 

-1.47 

44.52 
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